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Abstract 

The total cross section for the pp — > ppn reaction at energies close to thresh- 
old is calculated in a covariant one-boson-exchange model. The amplitudes for 
the elementary BN — > Ntt° processes are taken to be the sum of s, u and t pole 
terms. The main contributions to the primary production amplitude is due to 
a a meson exchange, which is strongly enhanced due to the t pole term. The 
NN and 7rN final state interactions are included coherently. The effects due to 
the 7rp interactions in pure isospin \ and | channels are sizable but when taken 
in the appropriate isospin combination almost cancel out. Both the scale and 
energy dependence of the cross section are perfectly reproduced. 
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In a recent study of the pp —* ppn reaction at energies close to threshold, 
it is found that the angular distributions of the outgoing particles are isotropic, 
in agreement with the assumption that the reaction proceeds as a 33 Po — ► 31 5oSo 
transition^], [|. Further more, the energy dependence of the total cross section, par- 
ticle spectra and angular distributions are well reproduced by taking into account 
the pp final state interactions (FSI) only The 7rN interactions seem to play a minor 
role, in marked difference with the similar pp —>■ ppr] reaction, where the 77N force is 
found to be essential for explaining the energy dependence of the cross section [|3]. |J . 
Several model calculations || 0, || of S-wave pion production, which are based on a 
single nucleon and a pion rescattering mechanism, under estimate the cross section 
by a factor of 3-5. Inspired by a study of /3-decay in nuclei, which indicates that 
the axial charge of the nuclear system is enhanced by heavy meson exchanges and, 
by using a simple operator form of the NN potential, Lee and Riska|| have shown 
that meson exchange currents could explain the scale of the cross section. Horowitz 
et al.[|IO| performed similar calculations based on explicit one-boson-exchange (OBE) 
model for the NN interaction as well as for the evaluation of meson exchange contri- 
butions. The evaluation of these contributions depends on the virtuality of the meson 
exchanged and on how the on mass-shell amplitudes are extended into the off mass- 
shell region. There are several approach based on field theoretical models which allow 
such extension to be made but the results are model dependent. In the traditional 
phenomenological treatment]!, ^, [5], U, the off mass-shell and on mass-shell ampli- 



tudes are of the same order of magnitude. Recently, several groups[ll, O, O 



have concluded that the off mass-shell 7rN rescattering amplitude is enhanced with 
respect to the on mass-shell amplitude. Hernandez and Oset[[Ll|] by applying current 
algebra and PC AC constraints argue that this enhancement may bring the calculated 
cross section for the pp — > ppir reaction into agreement with experiment. More de- 
tailed momentum-space calculations by Hanhart et al.Jl2| confirm this enhancement 
of the off mass-shell 7rN scattering amplitude but conclude that it is still too short to 
explain the scale of the cross section. Park et al. |13[ and Cohen et al. [[L4[] have applied 
a chiral perturbation theory, including chiral order and 1 Lagrangian terms. They 
have shown that the off mass-shell 7rN scattering amplitude is enhanced considerably 
but has an opposite sign with respect to the on mass-shell amplitude. Because of 
this difference in sign, the rescattering term and the Born term contributing to the 
pp — > ppir reaction interfere destructively, making the theoretical cross sections much 
smaller than experimental values. This enforces the importance of the heavy meson 
exchange contributions. Part of these contributions was included in Ref . [[14] . 



In the present work we propose a covariant OBE model based on a two-step mech- 
anism, where a virtual boson B (B=7r , o , r] , p , u ...) created on one of the incoming 
nucleons, is converted into a tt° meson on a the second via a BN — > tt°N process 
(see diagram of Fig. 1). Similar to tt° electroproduction amplitudes, the amplitudes 
for the elementary BN — > Nn processes are taken to be the sum of s, u and t pole 
terms. In fact only for the "effective" a meson, a t pole term does contribute, giving 
rise to a strong enhancement of the cross section. With this model we obtain perfect 
agreement with the empirical cross-section data. Though differing in many details, 



the present work agrees with Horowitz et al.|lQfl) that the a meson exchange plays 
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an important role in pion production near threshold. Our model however is a natu 
ral extension of the relativistic OBE model used to analyze NN scattering data [|15| 



where a single nucleon mechanism is not possible and all contributions correspond 
to completely connected diagrams. In the present formulation of the model, we may 
use any form of the BN —>■ tc°N amplitude so that comparison with other model 
calculations can be performed directly. 

We use the following Lagrangian interaction: 

L = ig 7rNN N-f 5 rNn + ig riNN N-f 5 Nri + g (TNN NNa + 

N [gpNNl 1 * + t^cr^q^rNp,+ 9uN N NYNu, + 9 &nn NtN5 , (1) 

with obvious notation. The coupling constants and the meson masses are taken from 
a fit to NN scattering data and are listed in table 5 of Ref.[T3|]. Using pseudovector 
couplings for the pseudoscalar mesons in the expression above would not alter the 
principal conclusions from the present work. 

Following Refs.[0, ffl we write the transition amplitude in the form, 



23 



M^T^ , (2) 



where M23 is the primary production amplitude which describes the transition from 
a two-nucleon state to a three-body state of two nucleon and a neutral pion, and T^P 
is a FSI correction factor taken to be the on mass-shell elastic scattering amplitude 
of the pp7i° — > ppir transition. The validity of this approximation for tc and rj meson 
production in NN collisions is discussed in length in Refs.|| [|. Here we note that : 
(1) all inelastic interactions ( and hence the coupling to other channels ) are included 
into M( m ) while T33 depends on elastic interactions among the reaction products 
only, (2) as in the coherent formalism for three-body processes |16|, |17|j , the different 
two-body interactions in the final state contribute coherently and (3) the transition 
amplitude has an overall phase factor identical to that of T33 as required by general 
quantum mechanics constraints ||16| . 



To calculate the transition amplitude we apply covariant perturbation theory tech- 
niques. We assume that the reaction is dominated by the mechanism depicted in Fig. 
1 and write M( m ' in the form, 

M (m) = £ [T BN ^ 7r0N (p 4 ,k;p 2 ,q)G B (q)S BNN (p 3 ,p 1 )] + [1 <-> 2; 3 <-> 4] , (3) 

B 

where pi, q and k are 4- momenta of the ith nucleon, the exchanged boson and the 
outgoing pion. The sum runs over all possible B boson exchanges that may contribute 
to the process. The bracket [1 <-» 2; 3 <-> 4] stands for a similar sum with the p±, p% 
and P2, Pa momenta interchanged. In Eqn. 3, T BN ^ n o N is the amplitude for the 
BN — > n°N transition, G B (q) and S B nn{P3,Pi) are the propagator and source func- 
tion of the meson exchanged, respectively. By using covariant meson propagators and 
form factors and covariant parametrizations for S B nn and T B ^^o^, the transition 
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amplitude is expressed in terms of invariant functions. For a scalar and pseudoscalar 
meson exchanges, 



Ssnn(Pi,P3) = u{p 3 ) I u(pi) F s (q) ; SW^pi, p 3 ) = "(Pa) 7 5 I w(pi) ^p(<?) , (4) 



where Fsiq) is a source form factor and / the appropriate isospin operator. Here u 
is a nucleon Dirac spinor and p% = pi — q is the final nucleon momentum. A most 
general vector source contains vector and tensor current terms, 



Svnn(PuPs) 



u( P3 ) YF^\q( 3 )+za^q u F^(q( 3 ) Iu{p 



(2), 



(5) 



with the vector source form factors, Fy , being analogous to the electromagnetic 
form factors of the nucleon. In the calculations to be presented below, all source form 
factors are taken in their lowest order perturbative approximation. The expression, 
Eqn. 5, contains conserved currents only and therefore satisfies current conservation. 

The amplitude for a scalar meson-nucleon scattering, SN — > PN, is written in 
the usual formflTSl, 



T. 



N(P4,k;p 2 ,q) = u(p 4 )r- 



U(P2 



(6) 



where A and B are functions of the Mandelstam variables, having the same isospin 
structure. To obtain the amplitude corresponding to a to a pseudoscalar meson 
scattering, PiN — > P 2 N, one has to remove the 7 5 in the expression above. 

Similarly, for a vector meson exchange , the amplitude for the VN — * tt°N tran- 
sition is, 



rpil 

± VN- 



■n°N 



(p 4 , k; j>2, q) = w(p 4 )7 5 [7 M ^i2 + PaMa + k"A m + q^A 7S \ u(p 2 



(7) 



where Aij stands for the combination (Ai + ftAj). 

In the analysis to be presented below, we include ir , a , r) , p and to meson 
exchanges. In order to calculate these contributions to the transition amplitude for 
the pp — > ppn reaction we call attention to certain kinematical features of this 
process. 

One important aspect of the mechanism of Fig. 1 is that at the ir° production 
threshold, the transferred 4-momentum is space-like, q 2 = — 3.3/m -2 . This is very 
much the same kinematics that occurs in ir° electroproduction through vector meson 
exchanges, suggesting that the half off mass-shell amplitudes Tbn^tt°n which appear 
in Eqn. 3, can be calculated using the usual calculation procedure of electroproduction 



amplitudes |]19||. Particularly, each of the amplitudes T v 



r)N^TT°N, T uN ^o N , lpN-*n°N 



and 



Tun-+k° n are far below the physical region and near the relevant poles, and therefore 
can be approximated as the sum of s, u and t pole terms. In fact, a t pole contributes 
to the To-jv-mOAT on ly an d as shown bellow it dominates the reaction amplitude. 

The situation is drastically different for 7r exchange because the T n o N^n°N is m the 
physical region close to threshold. It is now well established that the on mass-shell 
amplitude for the ir°p — > 7c°p scattering is very small and in the soft pion limit is equal 
to zeroPTJ near threshold. Based on general quantum mechanical constraints PO, 1TB 
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the half off mass-shell 7rN scattering amplitude can be written as the product of 
the on mass-shell amplitude and a real Kowalski-Noyes function|2(J, approaching 



unity in the limit of q 2 = m 2 . The derivation of this function is difficult and model 
dependent. In the traditional treatment of offshellness, this function corresponds to 
the form factor of the B meson leg (see diagram 1), which slightly enhances the off 
mass-shell amplitude with respect to the on mass-shell amplitude. Hernandez and 



Oset pl]] suggested two models that yield rather large enhancement. More detailed 



studies |12|, [H| [HJ] , agree that the enhancement due to offshellness, though differing 
in sign, is significant and must be taken into account. 

In order to study the influence of n exchange we evaluate the pp — > ppir cross 
section by (1) neglecting 71 exchange (this correspond to an exact cancellation between 
the direct and rescattering term of Refs.[TB|, fRJ]), and (2) using the enhanced 7rN 



amplitude of Park et al. . In agreement with Refs. [D| , it is found that, using 
the off mass-shell rescattering amplitude rather than the on mass-shell amplitude, 
influences the calculated cross section for the pp — > ppir rather little. 

As noted already, a significant contribution to the T a ^^ n o ^ amplitude is due to the 
t pole term. To calculate this require knowledge of the g anw coupling constant. Taking 
the Lagrangian density as L = g aiT7T M-n ■ na and by a simple one loop calculations 
this constant can be related to the g a NN through, 

ffL< ^ f^V djrNN. (\ n K _ {\ ^ ( ^NN^ ~~ 




4tt V 3 / 4:71 \ m n J \ An J 

where m n and M are the pion and proton mass. With the constants g^Aw/^ 71 " — 8-28 
and gl NN /4n = 14.6 one obtains g a7T7T = ( 3.3 ±0.1 ). The errors is due to neglecting 
terms of the order O y^kj- The transition amplitude of the reaction (solid curve) 
and the different meson exchange contributions are drawn in Fig. 2. The a meson 
exchange amplitude exceeds by far any other contribution. The relative phases for 
a , p , uj exchanges are predicted to be +1 and add constructively. The r\ exchange 
amplitude has a negative phase and interfere destructively, scaling down the transition 
amplitude to the proper magnitude required to explain the data. 
Finally, the cross section is calculated from the expression, 



(IT 



M 4 r d 3 p 3 d 3 p 4 d 3 p Ti 



16(2tt)V(s)Pi 



E 3 E4 E v 



Z \ 2 Sp (M (m) M (in)t ) 5\pi - Pf ), (9) 



where Sp (M( in ) M^t) denotes the trace over spinor states, and Z is the three- 
body FSI correction factor of Refs.|5], ffl. In the analysis presented below this factor 
is estimated from 7rN and NN elastic S-wave scattering phase shifts^, £|. The S-wave 
NN phase shift is obtained from the effective range expansion which includes Coulomb 
interaction between the two protons. We have used the scattering length a pp = —7.82 
fm and an effective range r pp = 2.7 fm of Ref.|23|. The Sll and S13 7rN scattering 



lengths are taken to be ai = 0.173 m^ 1 and 03 = —0.101 m' 1 , respectively |24| . 

Our predictions for the total cross section are drawn in Figs. 3 along with the 
available pp — > pp7T° data. The measured cross-section is reproduced remarkably well 
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with the 7r exchange contribution taken to be zero (solid curve). Note that all of the 
model parameters are determined independently from NN and 7rN scattering data 
and that none of these parameters have been adjusted to the pp — > ppn° reaction. 



Our predictions with the off mass-shell amplitude T^o p ^o p of Ref.|13| are drawn as 



a dashed curve. As the relative phase of the 7r exchange contribution and the other 
meson exchanges is not known we have repeated these calculations with the T n o p _> n o p 
taken with a reversed sign (small dashed line). Thus the effects of the 7r virtuality 
vary the calculated cross section by ~ 15% only, and as in Refs. ]T2|, [13|, [TJJ we may 
conclude that offshellness of the exchanged 7r would not resolve the discrepancy with 
experiment. 



In agreement with the model of Horowitz et 8l1. ||10|| , the main contribution is due 
to cr exchange with ratios M a : M w : M v : M p ps 100 : 40 : 8 : 7. Although the 
contribution of the r\ meson is rather weak it becomes effective through interference 
terms with the other exchange contributions. It is interesting to note that for the 
a exchange, the sum of the s and u pole terms alone amounts to ~ 6% of the r\ 
amplitude only, and that the main contribution is due to the t pole term which was 
not included in previous studies. The contribution from the 5 meson exchange is 
extremely weak for two reasons. The first is that the s and u pole terms have about 
the same magnitude but opposite signs. The second is that the mass of the 5 meson 
is high and the coupling constant gsNN is small. All these scale the 5 exchange 
contribution to less than 0.7% of the r\ exchange contribution. 

The transition amplitude is practically constant near threshold and the energy 
dependence is determined almost solely by the FSI correction factor and phase space. 
An important property of the coherent formalism is that the different two body in- 
teractions among the out going particles contribute coherently. Although the meson- 
nucleon interactions are weak with respect to the NN interaction, they become influ- 
ential through interference. To see this, we draw in Fig. 4 the cross section corrected 
for FSI interactions assuming pure l—h and I=| interactions for the outgoing 7r°N 
pairs. The energy dependence of the cross section differ significantly for the two 
channels. But, when the irN interactions are taken in the proper isospin combina- 
tion, their overall contribution to the FSI factor almost cancels out, leading to an 
energy dependence practically identical with the one corrected for the pp FSI (large 
dash curve) only. 

In summary, we have used OBE model to calculate S-wave production cross section 
for the pp — > ppn reaction by assuming a two-step mechanism where, a boson formed 
on one of the nucleons converts into a tt on the other through the BN — > 7r°N process. 
The kinematic conditions of the pp — ► ppn reaction at threshold, force the Tbn->tt°n 
amplitudes for the heavy mesons being far below the physical region, and hence could 
be estimated as the sum of s, u and t pole terms. It is found that the t pole term of o 
meson exchange is essential in order to reproduce the scale of the cross section. Using 
off mass-shell 7rN scattering amplitude as obtained from chiral perturbation theory 
would not affect these conclusions. 
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Figure 1: The primary production mechanism for the NN — > NNn reaction. 
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Figure 2: Predictions for the different meson exchange amplitudes to the pp — > ppir 
reaction vs. Q cm , the energy available in the center of mass system. The primary 
production amplitude is drawn as a solid curve. Note that the relative phases are 
taken into account. 
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Figure 3: Predictions for the total cross section of the pp — > ppir reaction. Predic- 
tions with the it exchange scattering amplitude taken to be zero are drawn as a solid 
curve. Those with the T v q p ^o p amplitude taken from Park et al. [[HJ are drawn as a 
dashed line. The small dashed curve gives predictions with the amplitude of Park et 
al.[P^] taken with the sign reversed. The data points are taken from Refs.|T], 
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Figure 4: FSI corrections in pure isospin 7rN channels. Integrated energy cross 
sections calculated with the assumption that the interacting 7rN pair is scattered in 
isospin I=| (small dashed curve) and I =| (dot-dashed curve). The solid line is that 
obtained with the ir N interactions taken in the appropriate isospin combinations. 
Predictions which account for the pp FSI only are (large dashed curve) are nearly 
identical with the solid curve. 
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